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Correlation between asymmetric profiles in slits and standard prewetting 

lines 

Salvador A. Sartarelli,-* Leszek Szybisd^'^^ 

The adsorption of Ar on substrates of Li is investigated within the framework of a den- 
sity functional theory which includes an effective pair potential recently proposed. This 
approach yields good results for the surface tension of the liquid- vapor interface over the 
entire range of temperatures, T, from the triple point, Tt, to the critical point, Tc. The be- 
havior of the adsorbate in the cases of a single planar wall and a slit geometry is analyzed as 
a function of temperature. Asymmetric density profiles are found for fluid confined in a slit 
built up of two identical planar walls leading to the spontaneous symmetry breaking (SSB) 
effect. We found that the asymmetric solutions occur even above the wetting temperature 
in a range of average densities Pg^^i < pav ^ p''ssb2; which diminishes with increasing 
temperatures until its disappearance at the critical prewetting point T^pw In this way a 
correlation between the disappearance of the SSB effect and the end of prewetting lines 
observed in the adsorption on a one-wall planar substrate is established. In addition, it 
is shown that a value for T^pw can be precisely determined by analyzing the asymmetry 
coefficients. 



I. Introduction 

The study of physisorption of fluids on solid sub- 
strates had led to very fascinating phenomena 
mainly determined by the relative strengths of 

*E-mail: asarta@ungs.edu. ar 

t E-mail: szybisz@tandar.cnea. gov. ar 

^ Institute de DesarroUo Humano, Universidad Nacional 
de General Sarmiento, Gutierrez 1150, RA-1663 San 
Miguel, Argentina. 

^ Laboratorio TANDAR, Departamento de Fi'sica, 
Comision Nacional de Energia Atomica, Av. del 
Libertador 8250, RA-1429 Buenos Aires, Argentina. 

^ Departamento de Fi'sica, Facultad de Ciencias Exactas y 
Naturales, Universidad de Buenos Aires, Ciudad Univer- 
sitaria, RA-1428 Buenos Aires, Argentina. 

* Consejo Nacional de Investigaciones Cienti'ficas y 
Tecnicas, Av. Rivadavia 1917, RA— 1033 Buenos Aires, 
Argentina. 



fluid-fluid (/-/) and substrate-fluid (s-/) attrac- 
tions. In the present work we shall refer to two 
of such features. One is the prewetting curve iden- 
tified in the study of fiuids adsorbed on planar sur- 
faces above the wetting temperature Tyj (see, e.g.. 
Pandit, Schick, and Wortis [l]) and the other is the 
occurrence of asymmetric profiles of fluids conflned 
in a slit of identical walls found by van Leeuwen 
and collaborators in molecular dynamics calcula- 
tions 2 , 3J . It is known that for a strong substrate 
(i.e., when the s-f attraction dominates over the 
/-/ one) the adsorbed film builds up continuously 
showing a complete wetting. In such a case, neither 
prewetting transitions nor spontaneous symmetry 
breaking (SSB) of the profiles are observed, both 
these phenomena appear for substrates of moder- 
ate strength. 

The prewetting has been widely analyzed for ad- 
sorption of quantum as well as classical fiuids. A 
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summary of experimental data and tlieoretical cal- 
culations for ''He may be found in Ref. [4]. Studies 
of other fluids are mentioned in Ref. [5]. These 
investigations indicated that prewetting is present 
in real systems such as ''He, H2, and inert gases 
adsorbed on alkali metals. 

On the other hand, after a recent work of Berim 
and Ruckenstein |6J there is a renewal of the inter- 
est in searching for the SSB effect in real systems. 
These authors utilized a density functional (DF) 
theory to study the confinement of Ar in a slit com- 
posed of two identical walls of CO2 and concluded 
that SSB occurs in a certain domain of tempera- 
tures. In a revised analysis of this case, reported in 
Ref. [7], we found that the conditions for the SSB 
were fulfilled because the authors of Ref. [6] had 
diminished the s-f attraction by locating an extra 
hard-wall repulsion. However, it was found that 
inert gases adsorbed on alkali metals exhibit SSB. 
Results for Ne confined by such substrates were re- 
cently reported |8|. 

The aim of the present investigation is to study 
the relation between the range of temperatures 
where the SSB occurs and the temperature depen- 
dence of the wetting properties. In this paper we 
illustrate our findings describing the results for Ar 
adsorbed on Li. Previous DF calculations of An- 
cilotto and Toigo [9] as well as Grand Canonical 
Monte Carlo (GCMC) simulations carried out by 
Curtarolo et al. \W suggest that Ar wets Li at a 
temperature significantly below Tc- So, this sys- 
tem should exhibit a large locus of the prewetting 
line and this feature makes it very convenient for 
our study as it was already communicated during 
a recent workshop [TT] . 

The paper is organized in the following way. The 
theoretical background is summarized in Sec. III. I 
The results, together with their analysis, are given 
in Sec. IIII.I Sec. llV.lis devoted to the conclusions. 



FDF[p{r)] 

j^id ksT Jdvp{v) {ln[A3p(r)] - 1} 

J rfrp(r)/Hs[p(r);dHs] 

J drdr/p{r)p{r//) $attr(| r - r/ |) 



+ / drp{r)Usfir) 



(1) 



The first term is the ideal gas free energy, 
where kg is the Boltzmann constant and A = 



2 TT h^/m ks T the de Broglie thermal wavelength 
of the molecule of mass m. Quantity is a pa- 
rameter introduced in Eq. (2) of [13] (in the stan- 
dard theory it is equal unity). The second term 
accounts for the repulsive /-/ interaction approxi- 
mated by a hard-sphere (HS) functional with a cer- 
tain choice for the HS diameter dns ■ In the present 
work we have used for /Hs[p(r); dns] the expression 
provided by the nonlocal DF (NLDF) formalism de- 
veloped by Kierlik and Rosinberg [M] (KR) , where 
p(r) is a properly averaged density. The third term 
is the attractive /-/ interactions treated in a mean 
field approximation (MFA). Finally, the last inte- 
gral represents the effect of the external potential 
Usf{r) exerted on the fluid. 

In the present work, for the analysis of ph- 
ysisorption we adopted the ab initio potential of 
Chismeshya, Cole, and Zaremba (CCZ) ^I5j with 
the parameters listed in Table 1 therein. 

i. Effective pair attraction 

The attractive part of the /-/ interaction was de- 
scribed by an effective pair interaction devised in 
Ref. [S] , where the separation of the Lennard- Jones 
(LJ) potential introduced by Weeks, Chandler and 
Andersen (WCA) [16j is adopted 



II. Theoretical background 

In a DF theory, the Helmholtz free energy Fdf [p(r)] 
of an inhomogeneous fluid embedded in an external 
potential Usf (r) is expressed as a functional of the 
local density p{r) (see, e.g., Ref. [T^) 



$WCA, 




(^) -M 



r < r„ 
r > r„ 



(2) 
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Here r„i — is the position of the LJ mini- 

mum. No cutoff for the pair potential was intro- 
duced. The well depth iff and the interaction 
size aff are considered as free parameters because 
the use of the bare values eff/ks = 119.76 K and 
<Jff = 3.405 A overestimates Tc. 

So, the complete DP formalism has three ad- 
justable parameters (namely, Vid, Sff, and (T//), 
which were determined by imposing that at l-v co- 
existence, the pressure as well as the chemical po- 
tential of the bulk I and v phases should be equal 
[i.e., P{pi) = P{pv) and p[pi) = p{pv)]- The pro- 
cedure is described in Ref. [5J. In practice, we set 
c^HS = ^ ff and imposed the coexistence data of pi , 
Ptj, and P{pi) = P{Pv) = ^0 for Ar quoted in Ta- 
ble X of Ref. 17J to be reproduced in the entire 
range of temperatures T between Tt — 83.78 K and 
Tc = 150.86 K. 



ii. Euler-Lagrange equation 

The equilibrium density profile p[r) of the adsorbed 
fluid is determined by a minimization of the free 
energy with respect to density variations with the 
constraint of a fixed number of particles N 



HSJ 



6piz) 

^ A , f A '^•^Hs[p(2/);rfHs] 5p(.zi) 

dzi p{zi) , , r , , . (6) 

'0 Sp{zi) 5p{z) 

Here Fid /A and Fhs/^ are free energies per unit of 
one wall area A. L is the size of the box adopted 
for solving the E-L equations. The boundary condi- 
tions for the one-wall and slit systems are different 
and will be given below. The final E-L equation 
may cast into the form 



i^id kB T In [A3 p{z)] + Q{z) = p 



(7) 



where 



Q{z) = f\is[p{z);d^s\ 

dzip{zi) — — — — 

5p{zi) 5p{z) 



+ 
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+ U,f{z). (8) 

The number of particles Ns per unit area, A, of the 
wall is 



Spir) 



FuF [p{r)] - p dv p{v) 



= 



(3) 



Here the Lagrange multiplier p is the chemical po- 
tential of the system. In the case of a planar sym- 
metry where the flat walls exhibit an infinite extent 
in the x and y directions, the profile depends only 
on the coordinate z perpendicular to the substrate. 
For this geometry, the variation of Eq. ([3]) yields 
the following Euler-Lagrange (E-L) equation 



<^[(Fid + Fhs)M] 
5p{z) 



where 
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i/id kB T In [k" p{z)] , (5) 



and 



N 



p{z) dz 
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In order to get solutions for p{z), it is useful to 
rewrite Eq. ([7|) as 



p{z) = pq exp 



Qjz) 
Vid kB T I 



with 



A-^ V i'id KB T 



(10) 



(11) 



The relation between p and Ns is obtained by sub- 
stituting Eq. pO)) into the constraint of Eq. ([9|) 



P = -Via kB T 
1 



X 



In 



iV.A3 



dz exp 



Qjz) \ 

i^idkBT J 



(12) 



When solving this kind of systems, it is usual to 
define dimensionless variables z* = zjbff for the 
distance and p* = p tr^j for the densities. In these 
units the box size becomes L* = Ljoff. 
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III. Results and Analysis 

In order to quantitatively study the adsorption of 
fluids within any theoretical approach, one must re- 
quire the experimental surface tension of the bulk 
liquid- vapor interface, 7;^, to be reproduced sat- 
isfactorily over the entire Tt < T < tempera- 
ture range. Therefore, we shall first examine the 
prediction for this observable before studying the 
adsorption phenomena. 

i. Surface tension of the bulk liquid-vapor 
interface 

Figure [T] shows the experimental data of 7;„ taken 
from Table II of Ref. [IF . In order to theoretically 
evaluate this quantity the E-L equations for free 
slabs of Ar, i.e. setting 

Usf{z) - , (13) 

were solved imposing periodic boundary conditions 
p{z = 0) = p{z = L). At a given temperature T, for 
a sufRciently large system one must obtain a wide 
central region with p{z ~ L/2) — pi{T) and tails 
with density pv{T), where the values of pi(T) and 
Pv{T) should be those of the liquid- vapor coexis- 
tence curve. The surface tension of the liquid- vapor 
interface is calculated according to the thermody- 
namic definition 

7i„ = (r! + Po V)/A = n/A + PoL, (14) 

where 51 = F^p — pN is the grand potential of the 
system and Pq the pressure at liquid-vapor coexis- 
tence previously introduced. We solved a box with 
L* ~ 40. The obtained results are plotted in Fig.[I] 
together with the prediction of the fluctuation the- 
ory of critical phenomena 7;^, = 7i°„(l — T/Tc)^'^^ 
with 7O = 17.4 K/A2 (see, e.g., [12] )■ One may 
realize that our values are in satisfactory agree- 
ment with experimental data and the renormaliza- 
tion theory over the entire range of temperatures 
Tt <T<Tc^ showing a small deviation near Tt. 

ii. Adsorption on one planar wall 

It is assumed that the physisorption of Ar on a one 
wall substrate of Li is driven by the CCZ potential, 
i.e.. 



12 




T[K] 

Figure 1: Surface tension of Ar as a function of 
temperature. Squares are experimental data taken 
from Table II of Ref. [H]. The solid curve corre- 
sponds to the fluctuation theory of critical phenom- 
ena and the circles are present DF results. 








Figure 2: Adsorption isotherms for the Ar/Li sys- 
tem, i.e., A/i as a function of coverage F^. Up- 
triangles correspond to T — 119 K; circles to 
T = 118 K; diamonds to T = 117 K; squares to 
T = 116 K; down-triangles to T = 114 K and stars 
to T = 112 K. 



Usf{z)^Uccz{z) . (15) 

The E-L equations were solved in a box of size 
L* — 40 by imposing p{z > L) = p{z = L). 
The solution gives a density profile p{z) and the 
corresponding chemical potential p. Adsorption 



010001-4 



Papers in Physics, vol. 1, art. 010001 (2009) / S. A. Sartarelli et al. 



isotherms at a given temperature were calculated as 
function of the excess surface density. This quan- 
tity, also termed coverage, is often expressed in 
nominal layers £ 

r£ = (i/p?^')/ dz[p{z)-pB], (16) 

Jo 

where pB — p{z ~> oo) is the asymptotic bulk den- 
sity and pi the liquid density at saturation for a 
given temperature. By utilizing the results for p 
obtained from the E-L equation and the value po 
corresponding to saturation at a given tempera- 
ture T, the difference Ap — p — pQ was evaluated. 
Figure [2] shows the adsorption isotherms for tem- 
peratures above T^, where an equal area Maxwell 
construction is feasible. This is just the prewet- 
ting region characterized by a jump in coverage Tg. 
The size of this jump depends on temperate. The 
largest jump occurs at and diminishes for in- 
creasing T until its disappearance at Tcpw Density 
profiles just below and above the coverage jump for 
T = 114 K are displayed in Fig. [3l in that case P^ 
jumps from 0.5 to 3.6. Therefore, the formation of 
the fourth layer may be observed in the plot. 




108 110 112 114 116 118 120 



T[K] 

Figure 4: Prewetting line for Ar adsorbed on Li. 
The solid curve is the fit to Eq. ([TT]) and reaches 
the Appw/ks = line at = 110.1 K. 

the jump in coverage occurs at each considered tem- 
perature. The behavior App^/kB vsT is displayed 
in Fig. m A useful form for determining the tem- 
perature Tu, was derived from thermodynamic ar- 
guments pp] 




App^{T) 



Figure 3: Examples of density profiles of Ar ad- 
sorbed on a surface of Li at T = 114 K displayed as 
a function of the distance from the wall located at 
z* ~ 0. Dashed curves are profiles for F^ below the 
coverage jump, while solid curves are stable films 
above this jump. 

The wetting temperature T^, can be obtained 
from the analysis of the values of Ap/kB at which 



fJ.pw{T) - po{T) 



(17) 



Here apw is a model parameter and the exponent 
3/2 is fixed by the power of the van der Walls tail 
of the adsorption potential Usf{z) ~ —C^/z^. The 
fit of the data of Ap/kB to Eq. ^ yielded = 
110.1 K and ap^/kB = -0.16 R-^/^. 

On the other hand, according to Fig. [21 the criti- 
cal prewetting point Tcpw lies between T = 118 and 
119 K. At the latter temperature, the film already 
presents a continuous growth. 

Our values of T^, and Tcpw are smaller than 
those obtained from prior DF calculations [5] {Tw = 
123 K and Tcpw ^ 130 K) and GCMC simulations 
[TD| (Tw = 130 K). The difference with the DF 
evaluation of Ref. [S] is due to the use of different 
effective pair potentials as we explain in Ref. [5], 
where the adsorption of Ne is studied. The present 
approach gives a reasonable 7i„ , while that of Ref. 
[5] fails dramatically close to Tf. The difference 
with the GCMC results cannot be interpreted in a 
straightforward way. 
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iii. Confinement in a planar slit 

In the slit geometry, where the Ar atoms are con- 
fined by two identical walls of Li the s-f potential 
becomes 

Usfiz)^Uccz{z) + UccziL-z) . (18) 

The walls were located at a distance L* =40, this 
width guarantees that the pair interaction between 
two atoms located at different walls is negligible. 
In fact, this width is wider than L* = 29.1, which 
was utilized in the pioneering molecular dynamics 
calculations [2] |3] . Accordingly, the E-L equations 
were solved in a box of size L* — 40. In this geom- 
etry, the repulsion at the walls causes the profiles 
p{z = 0) and p{z = L) to be equal to zero. The 
solutions were obtained at a fixed dimensionless av- 
erage density defined in terms of N ^ A, and L as 
pl, = Na)JAL = N*jL\ 



-11. 




Figure 5: Free energy per particle (in units of ks T) 
for Ar confined in a slit of Li with L* = 40 at 
T = 115 K displayed as a function of the average 
density. The curve labeled by circles corresponds to 
symmetric solutions, while that labeled by triangles 
corresponds to asymmetric ones. The SSB occurs 
in a certain range of average density /a*^^! ^ Plv ^ 

Pssb2- 

For temperatures below T^j = 110.1 K, we ob- 
tained large ranges of p*^ where the asymmetric 
solutions exhibit a lower free energy than the cor- 
responding symmetric ones. In spite of the fact 
that there is a general idea that a connection ex- 
ists between the SSB effect and nonwetting, we 



have found, by contrast, that SSB behavior extends 
above the wetting temperature. Furthermore, we 
have also found a relation between prewetting and 
SSB. 

Figure [5] shows the free energy per particle, 
/df = Fdf/N, for both symmetric and asymmetric 
solutions for the Ar/Li system atr=115K>Ttu 
as a function of the average density. According to 
this picture, the ground state (g-s) exhibits asym- 
metric profiles between a lower and an upper limit 
= 0.057 < pl, < = 0.192. Out of this 
range no asymmetric solutions were obtained form 
the set of Eqs. ((7)l- (fT2|) . Similar features were ob- 
tained for higher temperatures until T = 118 K, 
above this value the profiles corresponding to the 
g-s are always symmetric. Figure [5] shows three 
examples of solutions determined at T = 115 K. 
The result labeled 1 is a small asymmetric profile, 
that labeled 2 is the largest asymmetric solution at 
this temperature. So, by further increasing p*^, the 
SSB effect disappears and the g-s becomes symmet- 
ric, as indicated by the curve labeled 3. When the 
asymmetric profiles occur, the situation is denoted 
as partial (or one wall) wetting. The symmetric 
solutions account for a complete (two wall) wet- 
ting. These different situations can be interpreted 
in terms of the balance of 7^;, j^v a-nd 7;„ sur- 
face tensions, carefully discussed in previous works 
m |3l [7] . Here we shall restrict ourselves to briefiy 
outline the main features. When the liquid is ad- 
sorbed symmetrically like in the case of profile 3 
in Fig. [6l there are two s-l and two l-v interfaces. 
Hence, the total surface excess energy may be writ- 
ten as 

7,?r = 2 7.; + 2 7;.. (19) 
On the other hand, for a asymmetric profile j^^^ 
becomes 

itot = Isl + llv + Isv ■ (20) 

The three quantities of the r.h.s. of this equation 
are related by Young's law (see, e.g., Eq. (2.1) in 
Ref. [21]) 

Isv = Isl + llv COS.0 , (21) 

where 6 is the contact angle defined as the angle be- 
tween the wall and the interface between the liquid 
and the vapor (see Fig. 1 in Ref. [H]). By using 
Young's law, the Eq. (PD|) may be rewritten as 
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asy 
Itot 



2 7si + llv (1 + cos( 



(22) 



with cos6' = (7st, ~ ^si)hiv < 1- If one changes 
7sz by increasing enough Ng (as shown in Fig. [5]), 
and/or T, and/or the strength of Usf{z), eventually 
the equality jsv — Isi = liv may be reached yield- 
ing cos 9 — 1. Then, the system would undergo a 
transition to a symmetric profile where both walls 
of the slit are wet. 



1.4 



0.6 
0.4 
0.2 






Ar/Li T=115K 



10 15 20 25 30 35 40 



Figure 6: Density profiles of Ar confined in a slit 
of Li with L* = 40 at T = 115 K. The displayed 
spectra denoted by 1, 2 and 3 correspond to average 
densities p*„ = 0.074,0.192 and 0.218, respectively. 

It is important to remark that, indeed, there are 
two degenerate asymmetric solutions. Besides that 
one shown in Fig. [6] where the profiles exhibit the 
thicker film adsorbed on the left wall (left asym- 
metric solutions - LAS), there is an asymmetric so- 
lution with exactly the same free energy but where 
the thicker film is located near the right wall (right 
asymmetric solutions - RAS). 

The asymmetry of density profiles may be mea- 
sured by the quantity 



Af.s 



L/2 



dz [p{z) - p{L - z)] 



(23) 



According to this definition, if the profile is com- 
pletely asymmetrical about the middle of the slit, 
i.e. for: (i) p{z < L/2) ^ and p{z > L/2) = 0; 
or (ii) p{z < L/2) = and p{z > L/2) ^ this 




Figure 7: Asymmetry parameter for Ar confined 
by two Li walls separated by a distance of L* — 40 
as a function of average density. From outside to 
inside the curves correspond to temperatures T = 
112, 114, 115, 116, 117 and 118 K. The asymmetric 
solutions occur for different ranges p*gi,i < Pav < 

Pssb2- 




Figure 8: Circles stand for both branches of the 
asymmetry parameter for Ar confined in an L* = 
40 slit of Li walls for temperatures between and 



Tcpw The solid curve is the fit to Eq. (p4|) used to 
determine T, 



cpw • 



quantity becomes -1-1 or —1, respectively, while for 
symmetric solutions it vanishes. 

We evaluated the asymmetry coefficients of so- 
lutions obtained for increasing temperatures up to 
T = 118 K. The results for LAS profiles at tem- 
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peratures larger that are displayed in Fig. [7] as 
a function of the average density. One may ob- 
serve how the range < p*„ < p*,j2 dimin- 
ishes under increasing temperatures. The SSB ef- 
fect persists at most for the critical p*„(crii) = 
(17/24) CT^^ X 10~2 ~ 0.074 with ct// expressed in 

A. 

We shah demonstrate that by analyzing the data 
of A^r for p*^{crit) it is possible to determine the 
critical prewetting point. Figure [8] shows these 
values for both the LAS and RAS profiles, calcu- 
lated at different temperatures, suggesting a rather 
parabolic shape. So, we propose a fit to the follow- 
ing quartic polynomial 

T - r,p^ + a2A% + a4A% . (24) 

This procedure yielded Tcpw — 118.4 K, 02 = 
— 14.14 K, and a4 = —16.63 K. The obtained value 
of Tcpw is in agreement with the limits established 
when analyzing the adsorption isotherms of the 
one- wall systems displayed in Fig. [51 These results 
indicate that the disappearance of the SSB effect 
coincides with the end of the prewetting line. 

IV. Conclusions 

We have performed a consistent study within the 
same DF approach of free slabs of Ar, the adsorp- 
tion of these atoms on a single planar wall of Li 
and its confinement in slits of this alkali metal. 
Good results were obtained for the surface ten- 
sion of the liquid-vapor interface. The analysis of 
the physisorption on a planar surface indicates that 
Ar wets surfaces of Li in agreement with previous 
investigations. The isotherms for the adsorption 
on one planar wall exhibit a locus of prewetting 
in the /i — T plane. A fit of such data yielded a 
wetting temperature — 110.1 K. In addition, 
these isotherms also show that the critical prewet- 
ting point Tcpw lies between T — 118 and 119 K. 
These results for Tw and Tcpw are slightly below the 
values obtained in Refs. [9l[T0], the discrepancy is 
discussed in the text. 

On the other hand, this investigation shows that 
the profiles of Ar confined in a slit of Li present 
SSB. This effect occurs in a certain range of average 
densities p*^f,i < Pay < P*ssi,2J which diminishes for 
increasing temperatures. The main output of this 



work is the finding that above the wetting temper- 
ature the SSB occurs until Tcpw is reached. To the 
best of our knowledge this is the first time that such 
a correlation is reported. Furthermore, it is shown 
that by examining the evolution of the asymmetry 
coefficient one can precisely determine Tcpw- The 
obtained value Tcpw = 118.4 K lies in the inter- 
val established when analyzing the adsorption on a 
single wall. 
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